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A transient analysis for vesicle deformation under DC electric fields is developed. The theory extends from 
a droplet model, with the additional consideration of a lipid membrane separating two fluids of arbitrary 
properties. For the latter, both a membrane-charging and a membrane-mechanical model are supplied. The 
vesicle is assumed to remain spheroidal in shape for all times. The main result is an ODE governing the 
evolution of the vesicle aspect ratio. The effects of initial membrane tension and pulse length are examined. 
The model prediction is extensively compared with experimental data, and is shown to accurately capture the 
system behavior in the regime of no or weak electroporation. More importantly, the comparison reveals that 
vesicle relaxation obeys a universal behavior regardless of the means of deformation. The process is governed 
by a single timescale that is a function of the vesicle initial radius, the fluid viscosity, and the initial membrane 
tension. This universal scaling law can be used to calculate membrane properties from experimental data. 



I. INTRODUCTION 

Vesicles are widely used as a model system for biolog- 
ical cells due to their simplicity and controllability. The 
deformation of the lipid membrane, in particular under 
an applied electric field (electrodeformation) , is often ex- 
plored to probe membrane properties^ and to detect 
pathological changes in cells 4 

In the past decade, vesicle electrodeformation has be- 
come a significant subject of study, and earlier work 
can be divided into two categories. In the first cate- 
gory, an alternating-current (AC) field is applied, which 
often induces stationary and small deformations i^&L 
Correspondingly, an electrohydrodynamic theory in the 
small-deformation limit was developed to interpret the 
data trends^ In the second, under direct-current (DC) 
electric fields, vesicles usually exhibit large and tran- 
sient deformations due to the large field strengths com- 
monly applied^— Recently, using high-resolution, high- 
speed optical imaging Riske and DimovaS acquired a 
large amount of data capturing the complex deformation- 
relaxation behavior of the vesicles. Although some qual- 
itative and scaling arguments were presented^ the data 
was not fully interpreted due to the absence of a pre- 
dictive model. Meanwhile, one of us (HL) experimen- 
tally examined vesicles in the large-deformation regime 
with aspect ratios reaching tenjiii A large-deformation 
theory was also presented, which provided quantitative 
agreement with the data therein. However, the model 
was semi-empirical in that the hydrodynamic problem 
was not rigorously treated, but followed an empirical ap- 
proach by Hyuga and co-authors i 11 ' 12 In general, a rig- 
orous and transient analysis needs to be developed to 
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understand the complex deformation-relaxation behav- 
ior, and to provide insights on the underlying physical 
processes. 

In this work, we develop a transient analysis for vesi- 
cle electrodeformation. The theory is derived by extend- 
ing our previous work on a droplet model^ with the 
additional consideration of a lipid membrane separating 
two fluids of arbitrary properties. For the latter, both a 
membrane-charging and a membrane-mechanical model 
are supplied. Similar to the droplet model, the main re- 
sult is also an ordinary differential equation (ODE) gov- 
erning the evolution of the vesicle aspect ratio. The ef- 
fects of initial membrane tension and pulse length are 
examined. The model prediction is extensively com- 
pared with experimental data from Riske and Dimova 5 
and Sadik et al.^ and is shown to accurately capture 
the system behavior in the regime of no or weak elec- 
troporation. More importantly, the comparison reveals 
that vesicle relaxation obeys a universal behavior, and 
is governed by a single timescale that is a function of 
the vesicle initial radius, the fluid viscosity, and the ini- 
tial membrane tension. This behavior is regardless of the 
means of deformation, either via AC/DC electric field, or 
via mechanical stretching. This universal scaling law is a 
main contribution of the current work, and can be used to 
calculate membrane properties from experimental data. 



II. THEORY 

The problem configuration is shown in Fig. [TJ Under 
the influence of an applied electric field, charges of op- 
posite signs are allowed to accumulate on the two sides 
of the membrane, which induces vesicle deformation and 
electrohydrodynamic flows both inside and outside the 
vesicle. We assume that the vesicle remains spheroidal 
in shape throughout the process. All notations, as well 
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as the prolate spheroidal coordinate system follow those 
from Zhang et al.J£ The surface of the prolate spheriod 
is conveniently given as 



a 
c 



(1) 



Here c = \/a 2 — b 2 is chosen to be the semi-focal length 
of the spheroidal vesicle, and a and b are the major and 
minor semi-axis, respectively. For the derivation below, 
we further assume that the volume of the vesicle is con- 
served. We subsequently obtain 



ro(l -Z 2 y 



6 = r (l-£o 



(2) 



Therefore, the vesicle geometry is completely character- 
ized by a single parameter, £o> which evolves in time along 
with deformation. The critical idea of the current anal- 
ysis is to express all variables, e.g., the electric potential 
and the stream function in terms of £o- I n what follows, 
we introduce both an electrical and a mechanical model 
for the membrane. An ODE for £o is obtained by apply- 
ing the stress matching and kinematic conditions. 



A. The electrical problem 

The electric potentials both inside and outside the vesi- 
cle are described by the Laplace equations: 

V 2 2 = V 2 e = 0. (3) 

However, at the membrane the matching conditions are 
modified: 



<7e 

h* 



i) 



dt 



G,, 



i), at £ = &.(4) 



Here C m and G rn denote the membrane capacitance and 
conductance, respectively, is a metric coefficient of the 
prolate spheroidal coordinate system. This membrane- 
charging model is commonly adopted by many previ- 
ous research^ - — The displacement currents from the 
electrolytes are not included, which approximation is 
valid when the Maxwell- Wagner timescale, Tmw = (ej + 
2e e )/(<Ti + 2er e ), and the charge relaxation timescale, 
T cr = e/o") are small when compared with the membrane- 
charging time, T c h — r C m (l/ai + l/2a e ), and the de- 
formation time, Td — ii e /e e E^. However, the last two 
times are in general comparable with each other. The 
first term on the RHS of Eq. Q represents capacitive 
charging of the membrane, which includes the effect of 
membrane deformation. However, the contribution from 
this effect is usually small, and is neglected in the current 
analysis for simplicity. Equation Q can be consequently 
reduced to 

Ce d(f> e _ <Jj_dj)j_ _ C m c d(4> e - 

he d£ ~ he~~dt[ ~~h e di 



i), at £ = f .(5) 



(a) 

Eq 




rj = const 



FIG. 1. (a) A schematic of the problem configuration. The 
original radius of the vesicle is ro- The conductivity is denoted 
by a, the permittivity is denoted by e, the viscosity is denoted 
by and the subscripts i and e denote intravesicular and ex- 
travesicular, respectively. The strength of the applied electric 
field is Eq. (b) The prolate spheroidal coordinate system. 



Equation ([5]) can be further simplified by considering dif- 
ferent stages of charging. In the first stage, the trans- 
membrane potential (TMP), V m = ((pi — e )^ = ^ o , grows 
continuously in magnitude, but the membrane is not per- 
meabilized. Under this condition, G rn is near zero, and 
Eq. ([I]) becomes 



C m cd((j>e - 4>i) 



dt 



at£ = £ - (6) 



In the second stage, the maximum TMP reaches the 
critical threshold, V c , for electroporation to occurJ^— 
The membrane becomes permeable to ions, and G m in- 
creases significantly to limit further growth of the TMP. 
In general, the exact values of V m and G m depend on 
the detailed electroporation conditions and variables such 
as pore density and pore area.— The solution usually 
requires a complex numerical simulation which is be- 
yond the scope of the theoretical analysis pursued in 
this paper. However, a comprehensive model study by 
Li and Lii*i£ showed that the maximum TMP remained 
at the critical level in the presence of the pulse post- 
permeabilization. In this work, we adopt an approximate 
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model for this stage. We assume that once the maximum 
value of V m reaches V c , it no longer grows and "freezes" 
in time. In addition, the membrane is completely perme- 
abilized, and Eq. (|5|) is replaced by 



cr e d(p e 



V m = Vc, at £ = (7) 



Note that electroporation only occurs for sufficiently 
strong electric fields, and Eq. ([7]) is not needed for some 
of the cases studied below where V c is never reached. Far 
away from the vesicle surface, the electric field is uniform 



V0 e = E z, 



at £ — > oo. 



(8) 



We also require that (pi remains finite at £ = 1 . For initial 
condition, we solve Eqs. ([3]) and ((6]) with V m = 0. 



The general solution of the electric potentials for both 
the exterior and interior of the vesicle can be obtained 
following a similar procedure outlined in Zhang et al.i^ 



= E Q r Q [_A£ + aQi(0b, 



(pi = E o r o 0£r). 



(9) 



(10) 



Here, Qi(£) is a lst-degree Legendre polynomial of the 
second kind. A = c/ro is the dimcnsionless semi-focal 
length. The coefficients a and /3 are again obtained by 
applying the matching conditions. In the absence of elec- 
troporation, they are given as 



P + oyA 

Qi(£o)^r' 



(11) 



?i(£o)oy 



df3 
d^' 



(eo)Qi(eo)-Q?(go)(l-^) ^o | r 2 

dr nA 



f _ Qi(So) \ dA AQ' 1 , (e )gi(Co) 

40 QiteoMo Q'i 2 ^o) 



§=»• 



(12) 



«(0) 



A£o(e>V - 1) 



<2l(£o)£oOV - Ql(£o)' 



/3(0) = -A + a(0)Q x (£o) 



(13) 



Here ay = (J e /<Ji is the conductivity ratio. t± = r^Cm/di 
is a membrane-charging time. T2 = ro/ie/To is a charac- 
teristic flow timescale. To is the initial membrane ten- 
sion introduced below. The dimensionless time r defined 
as t = £/t2 has been used. Note that the definition of 
these times slightly deviates from those used in Zhang et 
al*£ due to the difference between droplet and vesicle. 
However, ti remains formally the same by replacing 7 in 
Zhang et al. 13 with To. 

After the maximum value of V m reaches the critical 
threshold, electroporation occurs, a and /3 are calculated 
by Eq. ([7]) which yields 



a= -V c /(E r )-\t Q (<j r -l) 

Ql(£o) ~ <2i(£o)£oOV 



(14) 

The expressions for the normal and tangential electro- 
static stresses are found in Zhang et al^ and not re- 
peated here. 



B. The hydrodynamic problem 

In the regime of low-Reynolds-number flow, the gov- 
erning equation for the hydrodynamic problem can be 



rewritten in terms of the stream function, tp, as 

E 4 ip = 0. (15) 

Here, the expression for the operator E 2 can be found in 
Dubash and Mestel 2 ^ and Bentenitis and Krause^ The 
stream function is related to the velocity components as 



1 dip 



1 dtp 
h v h e dr] 



(16) 



h v and hg are metric coefficients of the prolate spheroidal 
coordinate system. At the membrane, u and v represent 
the tangential and normal velocities, respectively, and 
they are required to be continuous 



Ui, 



at £ = Co- 



(17) 



In addition, we prescribe a kinematic condition relating 
the membrane displacement to the normal velocity, 

^ = e ,,) = ^^«% as) 



At the membrane, the stress matching condition is given 
as: 

||r-n|| =f mem . (19) 
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FIG. 2. The relative increase of the apparent area, A, as a 
function of membrane tension, Th, for different values of ini- 
tial membrane tension, Tq. The inset shows the linear regime 
for larger I\ values. 



Here f mem is the surface force density arising from the 
vesicle membrane. The tensor r includes contributions 
from both the hydrodynamic and electrostatic stresses: 



r = -pi + ^(Vv + Vv 5 



eEE - -e(E • E)I. (20) 



C. The membrane-mechanical model 

The surface force density at the vesicle membrane es- 
sentially consists of two parts&2£ 



(21) 



Here f K is the surface force density induced by bending 
resistance. f r = 2THn — V S T is the surface force density 
induced by the membrane tension. H is the mean curva- 
ture, and r is the local membrane tension. We can easily 



verify that f K is several orders of magnitude smaller than 
f r , and is therefore not included in the current analysis. 
The local membrane tension, F, is calculated by assum- 
ing an effective tension which is uniform over the entire 
membrane^ 2 ^ An increase of the homogeneous tension, 
Th, from the initial tension, Tq, leads to an increase in 
the apparent membrane area>i*22""— 



k B T Th . Th - r 
A = In 1 

8lTK Tq K a 



(22) 



Here A is the increase in the apparent membrane area 
relative to the initial spherical state, 



A 



-1. 



K a is the elastic stretching modulus, k is the bending 
rigidity. Equation (|22[) indicates that Tq, k, and K a are 
the important parameters in determining membrane ten- 
sion, k and K a are usually constants for a specific vesicle 
type, and their values are often readily obtained from 
previous workj 1 ' 32 ! 33 On the other hand, T is specific to 
an individual vesicle, and its value can not be directly 
determined from experimental measurements. The re- 
lation between A and Th for different choices of To is 
shown in Fig. [2] When A is small, the membrane area 
increases through the flattening of the undulations, and 
Th shows an exponential correlation with A. When A is 
sufficiently large, a linear behavior is observed instead, 
and the membrane area increase is mainly due to elastic 
stretching. Moreover, a larger Tq always leads to a larger 
Th for the same value of A. 



D. General solution 

A solution for vesicle electrodeformation can be ob- 
tained by solving the governing equations of both the 
electrical and hydrodynamic problems, with the help of 
the matching conditions. The solution strategy is iden- 
tical to that presented in Zhang et al.^ with only dif- 
ferences in the detailed matching conditions for both the 
electric field and the interfacial forces. For brevity, only 
the final governing equation for £o is presented here: 



dr F 



Q 



N 



Ca E 
A 2 



r\ ( (C \ i r\ ^r/22(£o) + /23(£o) Ffc 

yjV/2l(4o) + QT f J! \ , 7T : - — /24(4o) 

Mr/wUoJ + /l5l4oJ J-0 

(A - aQ'^o)) 2 + (A - aQ^o)/^) 2 - 2(3 2 /e r 



Qi 



Ca E r 



(A - ag 1 ($o))(A - aQi(G>)/£o) - P 2 hr 
I 



(24a) 
(24b) 
(24c) 



The functions /i4(^o) ; /i5(Co)j /2i(Co) — f2i{£,o), and F detailed expressions are found in the Appendix. e r = 
are the same as those used in Zhang et a/.,— and the 
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e e /ei is the permittivity ratio. The factors Qjv and Qt 
again arise from the effects of the tangential and normal 
stresses, respectively. Ccle = r e e EQ/T is the modified 
electric capillary number. In the absence of electropo- 
ration, the coefficients a and 8 are given in Eqs. (fTTj) 
and (|12|) . Once the electroporation occurs, Eq. fj 14|) is 
used instead. Similar to the droplet model, an examina- 
tion of the three terms in the numerator of Eq. (|24a[) 
reveals the contribution from the normal stress, tangen- 
tial stress, and membrane tension, respectively. The bal- 
ance between these three terms determines the equilib- 
rium vesicle shape. The above equations are solved until 
the end of the pulse, t = t p . 

In the context of vesicle electrodeformation, the relax- 
ation process is equally important, and is more reveal- 
ing of the underlying physical processes. The governing 
equations are presented below. In the absence of electro- 

I 



poration, Eq. ([3]) is solved without an applied electric 
field. The resulting equation for £o remains the same as 
Eq. (|24aj) . The coefficients of Qn, Qt, a, and B are 
given as 



IN 



A 2 r r 



a 2 (Q'i 2 (£o) + Q 2 (£o)/£o 2 ) -2/3 2 /e r 



A 2 r r 



(25) 



a 2 Qi(eo)Q 1 (^o)/eo-/3 2 /e I . , (26) 



Q'i(£o)oy ' 



(27) 



?i(£oK 



d3 
rf7 



l(£o)Qi(£o)-Q'i 2 (£o)(1 



dr 



T 2 
TlA 



(3 = 0, 



(28) 



a(T p ) 



V m (T p ) 



^o(Qi(6)&o-r-Qi 



f3(r P ) 



Wi(&)eoov-Oi 



(29) 



r 



In Eq. (|29p . the initial conditions for a and are ob- 
tained by solving Eqs. ([3]) and ([6]), and requiring that 
V m assumes the value at the end of the pulse. r p is the 
dimensionless time, t p /T2- Note that in this case, al- 
though the pulse is switched off, the electric field is in 
general not zero, due to the capacitive discharging of the 
membrane. In this case, the TMP will decreases from its 
peak value to zero on the membrane-charging timescale, 
T c h- 

When electroporation is present, the discharging pro- 
cess is slightly more complex. The full membrane- 
charging model (JSJ is used. In order to determine the 
membrane conductance, G m , we simply assume that it 
remains unchanged from the moment the pulse ceases, 
namely, 



(30) 



The resulting equation for £o again does not formally 
deviate from Eq. (|24a|) . The coefficients of Qn, Qt, ct, 
and are 



>N 



e e V c 2 
A 2 r r 



" 2 (Qi 2 (£o) + Q 2 (£o)/£ 2 ) -2/? 2 A 



QK 2 
A 2 r r 



(31) 



« 2 Qi(£o)Qi(£o)/£o-/3 2 Ar , (32) 





Q'i(£o)ov' 



(33) 



Qi(fr) 
Q'i(£o)ov 



ft 



d3_ 

dr 



Qi(eo)<9i(eo)-Qi 2 fe)(i-^ r )^o 



Q?(Zo)<Tr 



dr 



t 2 _ r 2 G m ( Qi(Co) 



13 = 0. 



(34) 



E. A similarity solution for vesicle relaxation we may ignore the membrane-discharging process. The 

membrane-charging/discharging time, T c h, is on the or- 
The governing equation for the relaxation process can der of 1 ms, which is in general much shorter than the 
be further simplified following two considerations. First, 
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relaxation time observed in the experiments, namely, a 
few tens of ms or longer. The relatively small effect of 
discharging on relaxation is clearly seen in Fig. |3] pre- 
sented in the following section. Without including the 
discharging process, the coefficients Qt and Qn in Eq. 
(|24a|) are simply set to zero. Second, in the membrane- 
mechanical model (|2"2"T) . the first and second term on 
the RHS represent the effects of undulation unfolding 
and elastic stretching, respectively. For moderate values 
of r , and for small-to- moderate deformations, the sec- 
ond term can be ignored, and the membrane-mechanical 
model becomes 



8ttk To 



Substituting Qt — Qn — and Eq. 
obtain 



d£ 1 8ttkA 



(35) 



33 into (I24al). we 



(36) 



This equation is conveniently rewritten in terms of the 
aspect ratio as 



ri- 
ll. 

dr 



1 ,8mA 2 
F^J^T^ 



(37) 



Note that in this equation, k, the bending rigidity, is 
regarded constant for a specific vesicle type, and /i r (em- 
bedded in F, see Appendix) is close to 1 as both the fluids 
are usually aqueous. In addition, A, the relative increase 
of apparent membrane area, depends exclusively on £oj 
hence a/b according to Eqs. (|23l) and ([2]). Under these 
assumptions, we observe that Eq. (1371) is completely au- 
tonomous, and the relaxation process is governed by the 
dimensionless time, t = t/r%, where Ti = ro/x e /Fo. This 
result suggests that the relaxation of vesicles with dif- 
ferent initial radius, r , and initial tension, r , obeys 
a similarity behavior with the proper scaling suggested 
above. This behavior is demonstrated by both simulation 
and analysis of previous experimental data below. 



III. RESULTS 

For all results below, we assume the lipid membrane 
to be made of egg-PC following Riske and Dimova 8 
(henceforth abbreviated as 'RD05') and Sadik et alJ£ 
(henceforth denoted as 'Sll'). The bending rigidity is 
taken to be k = 2.47 x 10~ 20 Jji the elastic modu- 
lus, K a = 0.14 N/m j 32 ' 33 the membrane capacitance, 
C m = 0.01 F/m 2 j21 the intravesicular and extravesicular 
viscosities, fa = fi e = 10 -3 Pa • s; the intravesicular and 
extravesicular permittivities, e, = e e = 7 x 10~ 10 F/m. 
The critical transmembrane potential is assumed to be 

V c = 1 VM- 



TABLE I. List of parameters for Fig. [5] For each case, Eq and 
t p are specified according to RD05. To is a fitting parameter 
to obtain best comparison between simulation and data. For 
cases b, d, e, and f, extended pulse lengths (denoted by star) 
are also used. 



case # 


E (kV/cm) 


tp (/is) 


To (N/m) 


a 


1 


150 


2.79 x 10" 4 


b 


1 


200 


3.23 x 10~ 6 




1 


300* 


3.23 x 10" 6 


c 


1 


250 


1.67 x 10" 4 


d 


1 


300 


1.80 x 10" 6 




1 


400* 


1.80 x 10~ 6 


e 


2 


50 


1.80 x 10" 4 




2 


80* 


1.80 x 10~ 4 


f 


2 


100 


3.16 x 10" 6 




2 


170* 


3.16 x 10~ 6 


g 


3 


50 


6.67 x 10 -6 


h 


3 


100 


3.42 x 10~ 7 



A. The effects of To and t p 

We begin by examining the effects of T on vesi- 
cle electrodeformation and relaxation. Figure [3] shows 
the typical system behavior for values of To ranging 
from 10~ 7 — 10~ 3 N/m. The intravesicular and ex- 
travesicular conductivities are Oi = 6 x 10 S/m and 
a e — 4.5 x 10~ 4 S/m, respectively following RD05. The 
field strength is Eq = 1 kV/cm, the pulse length is 
t p = 250 jtxs, and the initial radius is ro — 15 /im. Fig- 
ure EJa) shows the evolution of V m at the cathode- facing 
pole, which demonstrates only a weak dependence on Tq. 
The threshold for clectroporation (IV) is reached just 
before the end of the pulse, and its effects are present yet 
negligible. The discharging occurs on the relatively short 
timescale of 1 ms as we discussed above. Figure E^b) 
shows the evolution of the aspect ratio, a/b. The dis- 
charging process manifests itself as a sudden and slight 
decrease in the aspect ratio immediately after the pulse 
ceases; its effects can in general be ignored without signif- 
icantly altering the relaxation behavior. A smaller value 
of To leads to a larger aspect ratio, and a longer relax- 
ation process. The maximum aspect ratio, [a/6] max , is 
plotted as a function of To in Fig. [SJc). As the initial 
membrane tension decreases toward zero, the maximum 
achievable aspect ratio saturates. 

The similarity behavior in the relaxation process is 
demonstrated in Fig. ^d). The descending branches 
of the curves (t > t p ) shown in Fig. E^b) are rescaled in 
terms of r = t/r^, and shifted horizontally. In compari- 
son, the thick solid curve is obtained by directly solving 
Eq. (|3"T|) . The convergence of all curves validates that 
T 2 = roHe/To is the single timescale governing vesicle 
relaxation. 

The effects of t p are examined in Fig. 0] The pa- 
rameters are the same as in Fig. [3J and we fix Tq at 
1 x 10~ 6 N/m. Figure IDJa) shows that a longer pulse 
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FIG. 3. Vesicle deformation-relaxation as a function of To. The governing parameters are at — 6x 10 _4 S/m, a e = 4.5 x 10 _4 S/m, 
Eo = lkV/cm, t p — 250 /is, and ro = 15 /zm. (a) The transmembrane potential at the cathode-facing pole, (b) The time-course 
of the aspect ratio, (c) The maximum aspect ratio as a function of IV (d) The similarity behavior in relaxation. The descending 
branches from (b) are rescaled with r = t/T%. The thick solid curve is directly obtained by integrating Eq. (|37[) . 



consistently leads to greater deformation, and the aspect 
ratio increases along the same envelope. The relaxation 
times are approximately the same for all cases, because 
T2 remains unchanged. The discharging process is in gen- 
eral more conspicuous with longer pulses. In Fig. HJb), 
the relaxation curves are again shifted horizontally and 
rescaled with T2 to show good agreement with the simi- 
larity solution (thick solid line). Note that here because 
all cases share the same values of T2 , the collapse of the 
curves is primarily caused by simple shifting. In other 
words, the aspect ratio also decreases along a common 
envelope. 

The above results are exemplary and demonstrate the 
typical system behavior. In general, the relaxation pro- 
cess (in particular the relaxation time) is more apprecia- 
bly affected by the change in Tq than the deformation 



process. A wide range of pulsing parameters are studied 
below, in direct comparison with experimental data from 
RD05 and SIT. 



B. Comparison with experimental data 

An extensive comparison of our theoretical prediction 
with the data from RD05 is presented in Fig. [5j For 
all eight cases, the initial radius is ro = 15 fim. The 
electrical conductivities are cr, = 6 x 10~ 4 S/m and 
<j e = 4.5 x 1CP 4 S/m, respectively, leading to a conduc- 
tivity ratio of oy = 0.75. Other parameters are listed in 
table[U All parameters are taken directly from RD05, ex- 
cept for the extended pulse lengths for some cases noted 
below. For each case, the initial tension, Tq, is deter- 
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FIG. 4. Vesicle deformation- relaxation as a function of t p . 
The parameters are the same as in Fig. [3] The initial tension 
is set to be constant, To = 1 x 10 -6 N/m. (a) The time-course 
of the aspect ratio, (b) The similarity behavior is observed 
by shifting the relaxation curves with respect to time. The 
relaxation timescale, T2 = ro^i e /To, is the same for all cases. 
The thick solid curve is directly obtained by integrating Eq. 

03- 



mined to best fit the experimental data; their values are 
listed in table U in the last column. The experimental 
data are presented as symbols; the theoretical predic- 
tions, solid lines. In Figs. [S^a) to[S]Jd), the electric field 
strength is E = 1 kV/cm. For these cases, V m is pre- 
dicted to reach V c at t — 242 /is. In Figs. Eta) and 
[UJc) , good agreements are observed between the theoret- 
ical prediction and the data. In Figs. [5jb) andlSfd), the 
model results underpredict the maximum aspect ratios. 
This discrepancy is peculiar: our simulation follows the 
data accurately during the presence of the pulse, which 
duration is provided by RD05. After the pulse ceases, 



the simulation predicts immediate relaxation, whereas 
the vesicles continued to deform in the experiments, due 
to some unknown cause. In an attempt to mend this dif- 
ference, we artificially increase the pulse lengths in the 
simulation in b and d from 200 and 300 to 300 and 400 fis, 
respectively. The values for Tq remain unchanged. The 
results are shown as dashed curves. The model predicts 
well the data for both the deformation and relaxation 
processes. Note that although the relaxation curves rep- 
resented by the solid and dashed lines look somewhat 
different due to the semi-log scale on the time axis, they 
actually follow the same descending envelopes which we 
have demonstrated in Fig. QJb) above. 

In Figs. [5je) and Elf), the field strength is increased to 
be Eo — 2 kV/cm, and the pulse lengths used in RD05 
were 50 and 100 ^s, respectively. For these cases, our 
model predicts the occurrence of electroporation around 
t = 103 fj,s. A similar situation is observed as in Figs. 
Etb) and[5jd). The solid curves underpredict the maxi- 
mum aspect ratio. Artificially extending the pulses in e 
and f to 80 and 170 (is, respectively, leads to much better 
agreement between the two. 

In Figs. E{g) and[5jh), the field strength is further 
increased to 3 kV/cm, and electroporation is predicted 
to occur at t = 66 /is. The entire deformation-relaxation 
process is well-captured in g where t p = 50 /is. In Fig. 
E{h), where t p = 100 (is, although the model accurately 
predicts the deformation, the simulated relaxation curve 
completely deviates from the experimental data. For this 
case, and for pulses even longer than 100 /xs, RD05 [Fig. 
1(c) therein] exhibits a regime where complex, multi- 
stage relaxation process was observed. In this regime, 
the membrane structure is likely severely altered due to 
electroporation, which process can not be captured by 
our present model. Further comparison with these data 
is not pursued. 

The similarity behavior in the relaxation process is 
demonstrated in Fig. [6] The experimental data from 
Figs. Ela) to [S^g) are shifted horizontally and rescaled 
with T2- For each case, Ti is obtained using Tq listed in 
table |U The thick solid curve is again the similarity so- 
lution from Eq. (|37|) . and the results are shown on both 
semi-log and linear scales in r. The coefficient of de- 
termination is R 2 — 0.96. The experimental data from a 
wide range of parameters demonstrate a universal behav- 
ior governed by a single timescale, T2 = ro/ie/To- This 
result is a main contribution of the present work. 

We remark that a similar behavior should be observed 
for droplets, where the initial membrane tension, Tq, is 
replaced by 7, the coefficient of surface tension in t-i (cf. 
the definition of T2 in Zhang et al.). 13 However, there is 
a subtle difference between droplet and vesicle relaxation 
while the coefficient of surface tension is usually a con- 
stant, the membrane tension, I\, is not. Nonetheless, 
as long as Th depends linearly on To, which is a good 
approximation for small-to-modcratc deformations. The 
universal behavior in Fig. [6] is expected. 

Finally, the model prediction is compared with data 
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FIG. 5. Comparison with the deformation-relaxation data from RD05. For all cases, ro = 15 /im, Oi = 6 x 1(T 4 S/m, and 
<j e = 4.5 x 10~ 4 S/m. Parameters specific to each case are listed in table U The data is represented by symbols, and the 
simulation is represented by solid curves. For cases b, d, e, and f, the dashed lines represent the simulated results with extended 
pulses (denoted by stars in table HJ. 
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FIG. 6. The similarity behavior of vesicle relaxation. The 
experimental data from cases a-g in Fig. [5] are shifted in 
time, then rescaled by T2 = ro/ie/Fo. They are represented 
by symbols. The solid curves are calculated with Eq. (|37|) . 
The same data are shown on both a semi-log (a) and a linear 
(b) scale. The coefficient of determination is R 2 = 0.96. 



from Sll. In this work, the deformation is examined 
at a fixed pulse length of t p — 500 fjs, and for five 
intra-to-extra vesicular conductivity ratios. Only the 
case of Eq = 0.9 kV/cm is examined, where no or weak 
electroporation is expected. We do not compare the 
cases of Eq = 2 and 3 kV/cm in Sll, where the vesi- 
cles were in the strongly-electroporated regime, and our 
model no longer applies. The governing parameters are 
ro = 11.3 fim and a e — 3 x 10 -4 S/m. The initial mem- 
brane tension is chosen to be the same for all vesicles, 
namely, To = 1 x 10 -8 N/m. Figure[7][a) shows the defor- 
mation process as a function of time for five conductivity 
ratios. As ay decreases the rate of deformation increases. 
Except for the case of ay = 0.5, the aspect ratio reaches 



(a) 2.2 
2 
1.8 

ei-o 1 .6 
1.4 
1.2 



1 



-cr r = 5.00 X 10" 1 
■ a r = 1.67 x 10" 2 
a r = 3.13 x 10 -2 
cr r = 2.17 x 10" 2 
•a r = 1.89 x 10 -2 




100 200 300 400 500 

t (jjb) 



(b) 



2.5 



1.5 



Current Model 

* Sll 



20 



40 



60 



FIG. 7. Comparison with data from Sll. (a) Simulated time- 
course of the aspect ratio for various conductivity ratios. For 
all cases r Q = 11.3 ftm and r = 1 x 10~ 8 N/m. (b) The 
aspect ratio at t = 500 /is as a function of 1/oy. 



a plateau before the pulse ends. The time at which the 
aspect ratio increases saturates with an increasing ay. 
For ay = 0.5, an equilibrium could be reached if the 
pulse length is extended and sufficiently long (not shown 
here). In Fig. [TJb), the aspect ratio at t = t p is shown 
as a function of 1/oy. We choose this representation to 
facilitate comparison with the data from Sll (symbols), 
where the definition of the conductivity ratio is Ui/a e . 
A reasonable agreement is found between the two. The 
behavior of the simulation and the data is explained by 
the dependence of the electrical stress on ay in Sll [see 
Eq. (21) and Sec. 4 therein]. We do not repeat it here 
for brevity. The current model represents a significant 
improvement from that in Sll, where the hydrodynamic 
problem is treated empirically. 

Some remarks are appropriate before concluding the 
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section. First, for most cases studied here, the TMP is 
near the threshold, and the vesicles are expected to ex- 
perience no or weak electroporation. For this regime, 
our model is shown to provide a good predictive capabil- 
ity, which demonstrates that the membrane-mechanical 
model (|22|). although derived assuming no electropo- 
ration, can be extended to the weakly-electroporated 
regime, presumably due to the absence of major struc- 
tural alterations. Our model is not applicable to the 
strongly-electroporated regime. Second, the universal 
scaling law in relaxation observed in Figs. O HI and [5] is 
expected to hold regardless of the means of deformation, 
e.g., via AC/DC electric fields, or via mechanical stretch- 
ing. Equation <f37|) is applicable to a wide range of relax- 
ation phenomena beyond electrodeformation. Third, the 
current work suggests that an extensive parametric study 
on vesicle electrodeformation-relaxation experimentally, 
in particular in the sub-critical regime where electropora- 
tion is avoided, can provide the benefit to further validate 
our model understanding. A systematic approach can be 
possibly developed based on this work to map membrane 
properties. 



Riske and Dimova£ and Sadik et a/.j— and is shown to 
accurately capture the system behavior in the regime of 
no or weak electroporation. More importantly, the com- 
parison reveals that vesicle relaxation obeys a universal 
behavior, and is governed by a single timescale that is a 
function of the vesicle initial radius, the fluid viscosity, 
and the initial membrane tension. This behavior is re- 
gardless of the means of deformation, either via AC/DC 
electric field, or via mechanical stretching. This universal 
scaling law is a main contribution of the current work, 
and can be used to calculate membrane properties from 
experimental data. 
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Appendix A 



IV. CONCLUSIONS 

In this work, we developed a transient analysis for vesi- 
cle electrodeformation. The theory is derived by extend- 
ing our previous work on a droplet model in Zhang et 
al.r^ with the additional consideration of a lipid mem- 
brane separating two fluids of arbitrary properties. For 
the latter, both a membrane-charging and a membrane- 
mechanical model are supplied. Similar to the droplet 
model, the main result is also an ODE governing the 
evolution of the vesicle aspect ratio. The effects of initial 
membrane tension and pulse length are examined. The 
initial membrane tension affects the relaxation process 
much more significantly than the deformation process, 
in particular when its value is small. The model predic- 
tion is extensively compared with experimental data from 



The functions /u(Co), /is(Co), /2i(Co) - hi(£o), and 
F in Eq. (|24ap are given in the following expressions: 

/ute>)= / ^I drj, (Al) 
J (Co - T) 

(A2) 

t (c \ G 3 (Co)G 5 (Co) ~ G 3 (£p)G 5 (Co) , (c > . 
/l3lCo] = 7ii(,Co), (AdJ 

/i 4 (Co) = -Coff 3 (Co) J { ^ r ' ) r] L dr ] + ^ff 3 '(Co)/ii(Co), 

(A4) 



/is(Co) = - 
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G 3 {v)v 



drj. 



(A5) 



Here N = G 3 (Co)G 5 (Co) — G 3 (Co)G5(Co)- G and H are Gegenbauer functions of the first and second kind, respectively. 
The detailed expressions for G and H are found in Dassios et 
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M&H^l-^)* / L dv + Zo(l-Zo 2 T h I % = d V , (A9) 

i ? = -|(/26(&) + /26Ko)/Mr) ) (A10) 
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